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ABSTRACT 

We have examined the spatial distribution of substructure in clusters of galaxies us- 
QN ing Einstein X-ray observations. Subclusters are found to have a markedly anisotropic 
distribution which reflects the surrounding matter distribution on supercluster scales. 
! Our results suggest a picture in which cluster formation proceeds by mergers of sub- 
| £3 clusters along large-scale hlaments. The implications of such an anisotropic formation 
process for the shapes, orientations and kinematics of clusters are discussed briefly. 

<N 

Subject headings: cosmology: observations — galaxies: clusters: general — large-scale 
^ structure of the universe — X-rays: galaxies 

t> 
o 
m 

accepted for publication in the Astro-physical Journal Letters 

6 



1 Present address: Department of Astronomy &; Physics, Saint Mary's University, Halifax, NS B3H 3C3, Canada 
2 Sterrewacht Leiden, Postbus 9513, 2300 RA Leiden, The Netherlands 

3 Harvard- Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 



i 



1. Introduction 

Substructure appears to be a common feature of 
many, perhaps most, clusters of galaxies. Current es- 
timates suggest that at least 30% - 50% of rich clus- 
ters exhibit multiple concentrations in their galaxy or 
gas distribution (e.g. Geller k Beers 1982; Dressier 
k Shectman 1988; Jones k Forman 1992; Mohr, 
Fabricant k Geller 1993; Salvador-Sole, Sanroma k 
Gonzalez-Casado 1993; Bird 1994; Escalera et al. 
1994; Stern et al. 1995; see West 1994a for a recent 
review). 

Because dynamical evolution is expected to rapidly 
erase substructure, its prevalence in rich clusters to- 
day strongly suggests that we are currently in the 
epoch of cluster formation. Quantitative analysis of 
sub clustering may therefore provide information on 
cluster and galaxy formation, evolution of the intra- 
cluster medium and cosmology. For example, a num- 
ber of authors have suggested that the frequency of 
substructure in clusters today places a strong lower 
limit on the cosmological density parameter Slo (Gunn 
k Gott 1972; Richstone, Loeb k Turner 1992; Kauff- 
mann k White 1993; Lacey k Cole 1993). Similarly, 
the abundance and mass function of sub clusters may 
yield information about the primordial spectrum of 
density fluctuations (West, Oemler k Dekel 1988) as 
well as providing insights to galaxy formation (e.g., 
Beers k Geller 1983; Dressier 1984; Merritt 1985; 
Tremaine 1990). 

In this paper we endeavor to extract clues about 
the cluster formation process from the spatial distri- 
bution of sub clusters. Motivation for this study came 
from a number of intriguing coincidences that we had 
noticed between the distribution of subclusters and 
the surrounding matter distribution on much larger 
scales. An example is provided by the Coma cluster 
and its environs shown in Figure 1. Long considered 
to be the archetype of a rich, relaxed cluster of galax- 
ies, Coma is now known to possess a number of dis- 
tinct subclusters. The supergiant elliptical galaxies 
NGC 4889 and NGC 4874 reside in two subclusters 
in the core of Coma (Fitchett k Webster 1987; Mel- 
lier et al. 1988; Davis k Mushotzky 1993; Vikhlinin 
et al. 1994; Stern et al 1995) and ROSAT observa- 
tions have revealed a third large subcluster associated 
with NGC 4839, as well as several other smaller sub- 
clusters (Briel, Henry k Bohringer 1992; White et 
al. 1993). As Figure 1 reveals, the distribution of 
subclusters within the Coma cluster shows a rather 



striking alignment with the surrounding large-scale 
galaxy distribution, in particular the filamentary fea- 
ture which defines the Coma-Abell 1367 supercluster. 

A similar example is provided by Abell 426, which 
resides in the well-known Perseus-Pisces superclus- 
ter. R OSAT observations (Schwarz et al. 1992) show 
two X-ray peaks which lie right along the prominent 
supercluster ridge. Another member of the Perseus- 
Pisces supercluster, the poor cluster AWM 7 (Albert, 
White k Morgan 1977), also has substructure which 
shares the same orientation as the supercluster fila- 
ment (Stern et al. 1995). 

These and other examples suggest the intriguing 
possibility that the distribution of subclusters may 
be correlated with the surrounding matter distribu- 
tion on much larger scales. If true, this would be an 
important clue about the way in which galaxy clusters 
formed. However a few anecdotal cases like Coma or 
Perseus are not sufficient to establish the reality of 
this effect. Therefore, to examine this question more 
fully we undertook a statistical study of the distribu- 
tion of substructures in a large sample of clusters. 

2. The Cluster Sample 

X-ray observations have provided a wealth of in- 
formation on the structure of clusters of galaxies (see 
Forman k Jones 1990 and Jones k Forman 1992 for 
reviews). The largest substructure study to date is 
that of Jones and Forman (1992,1995) who assem- 
bled a sample of 366 clusters of galaxies with redshifts 
less than 0.2 observed with the Einstein satellite. 
Of these, 208 have adequate signal-to-noise to allow 
a reliable classification of their X-ray morphologies. 
Thirty-seven of these show clearly distinguished mul- 
tiple components, while an additional 56 are elliptical. 
The subclusters generally have projected separations 
corresponding to less than l/i -1 Mpc. It should be 
noted that other clusters in this sample, particularly 
the more distant ones, may also have substructure 
which was not detected owing to the arcminute spa- 
tial resolution of the Einstein Imaging Proportional 
Counter or to the superposition of substructures along 
the line of sight. For the 93 elliptical or multiple com- 
ponent clusters, we used either the positions of their 
subclusters or their ellipticity to define a position an- 
gle on the sky. 
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3. Linear Subcluster Configurations 

Visual inspection of the Einstein X-ray images 
gives a strong impression that when multiple sub clus- 
ters are present in a cluster they frequently have a 
collinear distribution. An example is shown in Figure 
2. In order to quantify this impression, we performed 
a simple statistical test using the seven clusters in our 
sample which have three distinct subcluster compo- 
nents (none of the clusters had more than three sub- 
clusters). Three subclusters define a triangle, whose 
shape can range from isosceles to a straight line. One 
can define a measure of the triangle shape as 

S= (D max - D mt )/D mm (1) 

where D m ; n ,D max and D; nt are the minimum, max- 
imum and intermediate angular separations between 
each of the three subcluster pairs. This ranges from 
S = for an isosceles triangle to S = 1 for a linear 
configuration. 

We computed this statistic for each of the three- 
component clusters in our sample; the results are pre- 
sented in Table 1. Despite the small sample size, the 
observed distribution of S values differs very signifi- 
cantly from that expected for random subcluster con- 
figurations. Monte Carlo simulations were performed 
of 10,000 random subcluster triplets, and these in- 
dicate that the distribution of S values in Table 1 
has a probability less than 1% of being consistent 
with randomly arranged subclusters. Hence there is 
a strong tendency for linear arrangements of subclus- 
ters in clusters. Presumably this is related to the 
linear shapes of many clusters of galaxies (Rood & 
Sastry 1971; Struble k Rood 1987). 

4. Subcluster Orientations with Respect to 
Large-Scale Structures 

One can test the idea that subclusters infall along 
preferred directions by comparing the orientation of 
the projected separation vector between each sub- 
cluster pair with the surrounding matter distribution 
on larger scales. Because most clusters are too dis- 
tant for the surrounding galaxy distribution to be 
seen in existing surveys, we used Abell clusters to 
map the surrounding large-scale structure. Although 
sparser tracers of the large-scale matter distribution 
than galaxies, rich clusters are known to delineate the 
same superclusters (e.g., Oort 1983; Bahcall 1992). 

For each of the 93 clusters in our sample, the 



projected position angle defined by its component 
subclusters or ellipticity, <f> ss , was compared with 
the projected position angle <f> cc from the cluster 
to each neighboring Abell cluster within a distance 
d < 10 /i -1 Mpc (position angles were computed us- 
ing standard spherical trigonometric relations). The 
difference between these two position angles defines 
an acute angle 9, 

= |^ss - </>cc| (2) 

which is a measure of the tendency for the subclus- 
ter and cluster pairs to be aligned with one another. 
Note that there is no ambiguity in assigning a value 
of <j) ss for clusters with three subclusters, owing to 
their strongly linear arrangements discussed in the 
previous section. Spatial separations between clus- 
ters were computed assuming a pure Hubble flow with 
Ho = 100 /i km s _1 and qo = 0.5. Cluster redshifts 
were taken primarily from the compilation by Pea- 
cock & West (1992) and augmented with a number of 
other recent measurements (Fetisova et al. 1993; Dal- 
ton et al. 1994; Lauer & Postman 1994; Quintana & 
Ramirez 1995). Figure 3 illustrates the method more 
clearly, along with the results. 

If the orientation of subcluster pairs is indepen- 
dent of the surrounding distribution of Abell clusters, 
then 9 should be uniformly distributed between 0°and 
90° . What one sees in Figure 3 is that for separa- 
tions less than 10 /i -1 Mpc the observed distribution 
is strongly skewed toward small values of 9, which in- 
dicates that subcluster pairs tend to share the same 
orientation as the surrounding large-scale cluster dis- 
tribution. Of the 93 clusters in our sample, 43 have 
one or more neighboring clusters within d = 10 /i -1 
Mpc. A Kolmogorov-Smirnov test confirms the statis- 
tical significance of these results; the probability that 
the distribution seen in Figure 3 could be consistent 
with a uniform distribution of 9 expected for random 
orientations of subcluster pairs is only ~ 0.6%. It is 
certainly significant that the distribution of subclus- 
ters appears to "know" about the surrounding distri- 
bution of Abell clusters on large scales. 

5. Discussion 

Our results provide clear evidence of a connection 
between the distribution of subclusters in galaxy clus- 
ters and the distribution of neighboring clusters on 
scales of ~ 10 h' 1 Mp c or more. 
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The most natural interpretation of these results 
is that cluster formation proceeds by the merging 
of subclusters along large-scale filamentary features 
in the matter distribution. This finding is not en- 
tirely unexpected, as N-body simulations have shown 
that filamentary infall may play an important role in 
some cluster formation models (e.g., West, Villumsen 
& Dekel 1991; Katz & White 1993; van Haarlem k 
van de Weygaert 1993; West 1994b, and references 
therein). 

Such an anisotropic formation process has impor- 
tant implications for the shapes, orientations and 
kinematics of clusters. In particular, this may explain 
the observed tendency for the major axes of Abell 
clusters to be aligned with their large-scale environs 
(e.g., Binggeli 1982; West 1989; Rhee, van Haarlem 
k Katgert 1992; Plionis 1994; West 1994b). Built up 
by a series of subcluster mergers which occur along 
preferred directions, clusters of galaxies will naturally 
develop orientations that reflect the surrounding fil- 
amentary pattern of superclustering. This formation 
scenario would also lead to strongly anisotropic clus- 
ter velocity dispersions. 

In conclusion, we have established a connection be- 
tween the internal structure of clusters of galaxies on 
sub-Mpc scales and the surrounding large-scale mat- 
ter distribution on supercluster scales. These results 
provide an important new clue about the genesis of 
galaxy clusters, suggesting that cluster formation pro- 
ceeds via the anisotropic merging of subclusters along 
filaments. 
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Table 1 



SUBCLUSTER CONFIGURATIONS 



Cluster 


S 


Abell 98 


0.998 


Abell 119 


0.773 


Abell 514 


0.404 


Abell 1477 


0.999 


Abell 1750 


0.998 


Abell 2384 


0.996 


Abell 4067 


0.988 
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Figure 1. Top panel: the large-scale galaxy distri- 
bution in the region of the Coma cluster. To high- 
light features in the galaxy distribution, symbol sizes 
are proportional to local galaxy density. Circles de- 
note Abell clusters with redshifts z < 0.03. Bottom 
panels: the distribution of subclusters in the Coma 
cluster. Note how the subclusters share the same ori- 
entation as the large-scale filament which defines the 
Coma-Abell 1367 supercluster. 

Figure 2. Einstein X-ray image of Abell 1750, show- 
ing the linear arrangement of its three component 
subclusters. 

Figure 3. Histogram of the distribution of 9 values 
for cluster separations d < 10 /i -1 Mpc, along with a 
schematic illustration of the method used. 



